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SUR QUELQUES MODÈLES STATISTIQUES 
SUGGÈRES 
PAR L'ÉTUDE DE L'EFFET DE SCINTILLATION 


par 


A. BLANC-LAPIERRE 


Professeur à la Faculté des Sciences d'Alger 


I. CARACTÉRISTIQUES EXPÉRIMENTALES DE L'EFFET DE SCIN- 
TILLATION 


On sait que les deux causes principales du bruit de fond dans un ampli- 
ficateur sont l'agitation thermique des électrons dans les résistances et la 
structure discontinue de l'électricité. La première cause donne naissance 
à ce qu'on appelle l'effet JOHNSON. La deuxième à l'effet de grenaille. 


a) Bruit d'agitation thermique : Soit une résistance R dans laquelle ne 
passe aucun courant. Îl existe cependant, entre ses deux extrémités A et 
B une différence de potentiel que l'on peut attribuer à une force électromo- 
trice fictive (placée en série avec une résis- 
tance de valeur R supposée sans bruit) et à 


laquelle on attribueraïit un spectre uniforme À 1 À 
de densité rie 
d[v]=4kRT (121) A 
R Nip 
(spectre réparti sur le 1/2 axe des fréquences a 
positives), k constante de BOLTZMANN, T Fo 
température absolue (Voir figure 1-1). 3 ; à 


b) Effet de grenaille : Si on considère un ete, 
courant électronique 1 émis par une cathode Mu e(t) 
saturée, I se comporte comme une fonction 
aléatoire I(t) à laquelle serait associé un spec- 
tre de densité constante Fig. (1-1) 


d[v]=2el (1-2 
(1 = valeur moyenne = intensité au sens macroscopique, e = charge de 
l'électron) (spectre réparti sur le 1/2 axe des fréquences positives). 


Les deux effets précédents correspondent à des fluctuations de spec- 
tre uniforme (tout au moins pour toutes les fréquences utilisées). 
Or, l'expérience a montré que, dans certains cas, aux phénomènes de 


fluctuations précédents, se superposait un autre phénomène de fluctuation, 
intervenant seulement pratiquement dans le domaine des basses fréquences 


__— 
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et qui est lié au passage d'un courant permanent. Ce phénomène se présen- 
te dans des circonstances assez diverses mais ilse présente toujours avec 
des caractéristiques statistiques semblables. C'est cet accroissement de 
bruit en basses fréquences que, pour des raisons qui seront claires plus 
loin, nous appellerons effet de scintillation. 


Avant d'aller plus loin, donnons une idée des circonstances dans les- 
quelles on le rencontre : 


a) Considérons une résistance métallique. En l'absence de courant 
permanent elle présente des fluctuations de potentiel conformes à (1-1); 
faisons passer un courant permanent : s'il s'agit d'une résistance consti- 
tuée de métal massif, il n'y a rien de changé. Si, au contraire, il s'agit 
d'une résistance de carbone aggloméré [1] ou d'un métal en couche mince 
[2]1[3] ou de couches minces de carbone [4], il y a un accroissement du 
bruit en basse fréquence : c'est l'effet de scintillation. De même s'il s'agit 
d'un cristal détecteur [5]. 


b) Considérons une diode pourvue d'une cathode à oxydes [6]; si on 
opère dans un domaine de fréquences assez élevées, les fluctuations de 
courant sont conformes à (1-2). En basses fréquences (disons pour 
V £ 20 000 cs/S), on observe un bruit qui peut être bien supérieur à celui 
qui est indiqué par la formule (1-2). Cet effet est très marqué pour les ca- 
thodes à oxydes; on l'a, cependant, aussi observé sur des cathodes consti- 
tuées par des filaments de tungstène et la première interprétation qui en a 
été donnée a été la suivante : il s'agit de modifications locales et acciden- 
telles du pouvoir émissif de la cathode sous l'effet d'absorption d'atomes à 
la surface; la surface de la cathode, frappée par les atomes du gaz rési- 
duel, fixe certains de ceux-ci puis les libère au bout d'un temps aléatoire. 
C'est cette image qui a suggéré la dénomination d'effet de scintillation. 


Quelles sont les propriétés statistiques de l'effet de scintillation? 


Dans tous les cas, le phénomène obéit à certaines lois qui paraissent 
avoir une portée tout à fait générale. La densité spectrale est une fonction 
décroissante de la fréquence V (tout au moins si on excepte la zone relati- 
vement peu étudiée des très basses fréquences; par exemple, V < 50 c/s) 
qui dans le domaine où l'effet est important est du type 


d [v]=rvA4/v (1-3) 


- dans de très larges limites de fréquence, d[v] décroit avec la fré- 
quence suivant une loi du type (1-3) avec 


X # 1 


- pour les fréquences plus élevées, d {v] décroit plus vite; « peut alors 
devenir voisin de 2. D'ailleurs il faut bien assurer la convergence de 


J dv 


- par ailleurs d (v)est proportionnel à I" avec n très souvent voisin de 
2 - c'est ce que nous appellerons la "oi en I". 


7 Dans le cas des lames minces ou des cristaux détecteurs, N. NIFON- 
OFF a montré de façon expérimentale qu'une liaison très étroite existait 
entre l'écart à la loi d'Ohm et l'écart à la loi en I’. 


Tels sont, dans leur ensemble, les faits expérimentaux. 
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L'explication physique détaillée n'est probablement pas encore élabo- 
rée de façon définitive. D'ailleurs, selon qu'il s'agit de résistance de 1la- 
mes minces, de cristaux détecteurs ou de cathode à oxyde, le processus 
physique n'est pas le même. Par exemple, dans le cas des lames minces, 
il y a probablement fluctuation du nombre d'électrons libres dans chaque 
granule et fluctuation des barrières de potentiel entre les granules... Ce- 
pendant, le fait que l'allure statistique de tous ces effets est la même, 
conduit à penser que, indépendamment de toute interprétation physique 
trop précise, il est intéressant de développer systématiquement l'étude de 
modèles statistiques présentant des régimes de fluctuations ayant des pro- 
priétés voisines de celles de l'effet de scintillation. Au demeurant, de tels 
modèles sont probablement utiles dans de nombreuses questions extérieu- 
res à l'étude de l'effet de scintillation [7]. 


Avec des complications et des variantes diverses, tous les modèles 
que je vais décrire seront plus ou moins rattachés à des distributions de 
POISSON de sorte que l'on peut dire qu'ils constituent des exemples pris 
parmi la riche famille des fonctions aélatoires dérivées des distributions 
de POISSON. 


II. SCHÉMA N° 1 


Ce schéma avait été proposé par SCHOTTKY (8 ) à la suite des expé- 
riences de JOHNSON (6). SCHOTTKY en avait ébauché l'étude. 


Nous partons d'une distribution de POISSON de densité uniforme P,; tj 
sera un point de cette distribution. À chaque instant tj, nous associons une 
impulsion rectangulaire de valeur 1 et de longueur aléatoire tj t'j; nous 
supposons que les longueurs des impulsions sont des variables aléatoires 
indépendantes entre elles et que la probabilité pour qu'une impulsion ait 
une longueur 1 (à dl près) soit égale à 


e? dl (2-1) 


x est la durée moyenne des impulsions. La figure (2-1) rendra plus clair 
le mécanisme aléatoire décrit. Il n'y a aucune impossibilité à ce que les 
impulsions empiètent les unes sur les autres. Soit x(t) la somme des va- 
leurs correspondant à l'ensemble des impulsions non nulles à l'instant t. 
C'est x(t) que nous associons au processus de POISSON. 


Fig. (2-1) 


1 + 
t; tj bat tj+2 Ciss Lire 


L'interprétation physique de ce schéma est la suivante : considérons, 
par exemple, une cathode de tungstène, les instants tj seront ceux où une 
molécule (ou un atome) vient se fixer par adsorption sur la surface émis- 
sive. -L'instant t'j associé à tj est celui où la molécule fixée en tj quitte la 
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surface. À est la durée moyenne de séjour d'une molécule à la surface de 
la cathode et x représente le nombre de molécules présentées, à la surface, 
à l'instant t. SCHOTTKY admet que le courant d'émission de la cathode dif- 
fère de la valeur I qu'il aurait en l'absence de toute molécule adsorbée 
d'une valeur AI proportionnelle à I et à x(t). 


AI KI xt) (2-2) 
soit 
ERA} LES RAI TExE (2-3) 


J'ai donné [9] une étude directe de la fonction aléatoire x(t). Je pense 
qu'il est bon de rattacher l'étude de x(t) à celle d'une classe très générale 
de fonctions aléatoires; la classe des fonctions aléatoires qui dérivent d'une 
distribution de POISSON par une formule du type 


X-(E)2= 72 Et, 4) (2-4) 
où, pour tj donné, R (t, tj) est une fonction aléatoire de t, les termes 
R (t, tj) associés aux différents tj étant indépendants les uns des autres. 


R. FORTET a étendu à ces fonctions la méthode que j'avais utilisée pour 
l'étude de l'effet de grenaille pur. Il est parvenu au résultat suivant [10] : 


Soit [u,, ...u,;t,, t,...t,; t;] la fonction caractéristique relative aux 
k variables aléatoires : 


REC an er nl 
Posons d'autre part 
de CO) LS EX Dee OEM SX RTE (2-5) 


Alors la fonction caractéristique ® (u,... u,) relative à x,...x, est 
telle que : ro 


e- P LP Chosrsatrs 0) oil d0 (2-6) 


avec Y = log®. 


Limitons-nous dans l'étude du schéma actuel au cas k = 2 et calculons 


nue t., 0 ). Nous pouvons supposer t, < t,. Tenons alors compte de ce 


a) la probabilité pour qu'une impulsion ait long i 
PA pret q p it une longueur comprise en- 


Pæ=1-6"1h (2-7) 


b) la probabilité pour que cette impulsion ait une 1 
entre 1, et 1, est 1 P une longueur comprise 


Base AL en (1222079 (2-8) 
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c) la probabilité pour qu 


l'est 


On obtient alors, 


à t'et t: 
ANIME TE: 
(a — fre tra 
“ e MUsHiUs 


p = 1-e-ttr-0)/x 
Chant: SD 
PI 


D'où la fonction y 


P=rer ra 


lune impulsion ait une longueur supérieure à 


(2-9) 


pour , Suivant la façon dont 0 se situe par rapport 


Fos [ e“t1- 01/2 


e-(t:-0)J2 


ei" e {te -60)/2 


cherchée 


Ne =(t3 = 01) D 


(2-10) 


(2-11) 


y pe x ne + e'" ( Te EPS Le EU ea] 


- - i -(t2-t,)/X 
PAR D PET Cnil I 


(2-12) 


y — Po À fre DE Le LP NET) ie PR 


Si l'on s'intéresse non plus aux x mais aux 


et si on utilise 


ES le tro) 


X 


(Po 


xX'= 


on obtient, pour les grandes valeurs dep, , une loi de LA PLACE tangente 
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dont la fonction caractéristique est donnée par 


W'= Log &'[U,.U,] =-4 {0 +Ui + 2U, U, NUE (2-13) 


De plus, quelle que soit la valeur de Vp,\ , on a, pour la fonction de 
corrélation 


GTS ONE OUEN ER (2-14) 


Le spectre associé à cette fonction de corrélation est défini par (v) 
tel que 


GS) NL SEX (2-15) 
dv HHATTAVEO 


f (v) = 


Il admet une densité f (v). La figure (2-2) représente la fonction de cor- 
rélation et la densité dont nous venons de donner les expressions analyti- 
ques. Pour des valeurs de v assez grandes, onaune densité f(v) qui décroit 
comme 1/V°. 


À À 
C(G) | F(0) 


Fig. (2-2) 


Remarques - Avant de quitter l'étude de ce schéma il est bon de faire 
quelques remarques. 


1°) - Supposons que la densité de la loi de POISSON ne soit pas cons- 
tante mais, par exemple, périodique. Alors la fonction aléatoire X(t) n'est 
plus stationnaire. On obtient sa covariance E[X(t,).X(t,)] en cherchant le 
nie en U, U, dans le développement de U,U,; il provient uniquement du 
erme en 


iU, + iUz 


e 


de (2-10) 
On trouve 


ts 
ER UNE Î p(O)etr02 40 
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D'où 


t 
COUR - | p (0) et gp 


ei eo p (0) e® dO 


On sait que, dans ce cas, on peut encore définir un spectre à condition 
d'introduire la moyenne temporelle ( 11) (12) : 


lim nn EMULE EC(E +h9 Lt AT th) 


T+ oo 


qui jouera le rôle de la fonction de corrélation; or, on voit sans peine que 


+T + 
_lim + dt{ f Pope MAO Le Sn (2-16) 


où P estla moyenne temporelle dep. C'est en effet la valeur de la com- 
posante continue de la transformée de p (t) dans un filtre de réponse per- 
cussionnelle R telle que 


R=0 pour @G@<0 
R=e%* pour 0>0 
On a donc 
PCR xNe NI /2 (2-17) 


On retrouve la formule (2-15). Le spectre lié à la moyenne temporelle 
de la puissance est donc encore égal à celui qui est donné par (2-16). 


2°) - Bien que ceci nous écarte de l'effet de scintillation, faisons re- 
marquer que la classe des fonctions aléatoires (2-4) contient, comme cas 
particulier, des fonctions aléatoires introduites par l'étude de l'émission 


= 


secondaire. Reportons-nous à la figure (2-3). I, représente un courant 


Is 


Collecteur 


Paguets d ‘electrons 
secondaires 


Flectrons 


primaires X(E) 


( 


Cible 
à émission Secondaire Fig. (2-3) 


Amplificateur 
a) 
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électronique jouant le rôle de courant primaire. Ce courant primaire tom- 
be sur la cible C dont on utilise l'émission secondaire. Chaque électron e; 
du courant primaire donne lieu à l'apparition simultanée de n; électrons se- 
condaires. Les n; sont des variables aléatoires indépendantes. Le courant 
secondaire est donc constitué par des paquets de n; électrons. Si ce cou- 
rant secondaire est appliqué à l'entrée d'un amplificateur de réponse per- 
cussionnelle r, on observera à la sortie de ce dernier 


X (T)=S RE De Envie) 


R est aléatoire par n;. Si w' est la fonction caractéristique relative à n;, 
on déduit de (2-7) l'expression suivante qui permet l'étude de la plupart des 
phénomènes de fluctuations dans l'émission secondaire 


vf" [ete r(ti -6)}-1]d0 


HI. SCHÉMA N° 2 


Nous partons encore d'une distribution de POISSON commençant à 
t =-00, Ily aune probabilité nulle pour qu'il y ait un point situé exactement 
en t = 0. On peut donc parler sans ambiguité de l'intervalle tj, t;,, conte- 
nant l'origine. Tirons alors au sort entre deux valeurs R, et R, supposées 
équiprobables. Pour tj<t < tj:1, nous donnons à x(t) la valeur fournie par 
le tirage et nous convenons quià chaque instant t, correspondant à un saut 
dans la distribution de POISSON, on passe alternativement de l'une à l'au- 
tre des valeurs R, et R, (voir figure 3-1). Ce schéma peut avoir quelques 
rapports avec le problème suivant qui nous intéresse : s'il s'agit d'une ré- 


l 


Fig. (3-1) 
Ra 


Ri 


1 
| 
l 
1 
! 
l 
! 
1 
i 
! 
l 
! 
1 
[ 


4 EE: 1 mir Se 


je 
b Liu lj+2 b+s (4 


sistance constituée par du carbone aggloméré. En présence d'un courant 
permanent I nous pou7ons admettre que les résistances de contact subis- 
sent des modificaticns accidentelles et désordonnées. À une variation AR 
correspond une variation 


AV=IAR (3-1) 


de la différence de potentiel; le schéma actuel est un modèle très simplifié 
où R ne peut prendre que 2 valeurs. 
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Pour étudier les propriétés de x(t) nous supposerons que l'on a 


Ri=-1! et R,= +1 


1°) Fonctions de répartition - Fonctions caractéristiques 


.. Nous appelons dn (0) le nombre d'instants présentés par la distribu- 
tion de POISSON entre 0 et 0 + d0; le nombre moyen à de points tj tom- 
bant entre tt, et t, est 


œ=E{ f'an(0)}-f"p (0) à 0 (3-2) 


(pour l'instant, on ne suppose pas p constant). 


Cela posé, nous allons calculer la loi de probabilité qui régit l'ensem- 
ble des deux variables aléatoires liées 


x(t;:)=x, X (t,Ji=, x (3-3) 
La probabilité pour qu'il y ait un nombre pair de chocs entre t, ett, est 
P+ = e° ch (3-4) 
Celle pour qu'il y ait un nombre impair de chocs est 

P-= e * sh x (3-5) 
Si on considère seulement x(t), on a une probabilité 1/2 d'avoir 
X,=+1 et une probabilité 1/2 d'avoir x, = - 1. Ces remarques simples 
étant faites, il est facile de calculer la fonction caractéristique associée à 
l'ensemble (x,, x,). Pour calculer E ele il suffit de se reporter au 

tableau ci-dessous qui résume l'ensemble des cas possibles. On trouve 


[uv A = e”"*[ ch œ cos (U, +u,) + sh œ cos(u;-U:)] (3-6) 


| 


ca 
Prob T2 Pi 


EN. 


RH Prob =} Pi 


2°) Propriétés du second ordre 
On déduit de l'expression de la fonction caractéristique, l'expression 
de la covariance 


E[x(t,) x(t, Je e-A*Pto) de (3-7) 
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A - Cas d'une distribution de POISSON stationnaire 
Alors 


p (0) = p. = constante 
CHU) SEE) hs (EE TD 4] Ro El (3-8) 


Le spectre correspond à une densité 


FN) EEE (3-9) 


© 
PE en Aie VE 


qui est de forme analogue à la densité définie par (2-16) 


B - Cas d'une distribution de POISSON non stationnaire 


On peut encore avoir un spectre en introduisant les moyennes tempo- 
relles si elles existent. On fait alors jouer le rôle de fonction de corréla- 
tion à 


Déte le din Et mens (3-10) 


T—+ 00 2 ST 


On peut pousser les calculs jusqu'au bout sans difficulté si 


p (9)= Ps + p, cos. w,:0 (3-11) 
avec naturellement 
Pr © De 
On trouve 
F'érYsterentre Her ur (3-12) 


Si p, 0, I, +1 et on retrouve la relation (3-8) valable dans le cas station- 
naire. Si on développe la fonction 


AOL TT 
LA sin =) 
en série de FOURIER, on peut écrire 


To) _ e-2PIT B, eikwe T (3-13) 


avec 


M ONE D MU p ee (ar)! (3-14) 
‘ tn) UCI TE KE CE HIKID 
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La densité f(v) Se présentera sous la forme d'une somme de fonctions f,(v) 
définies par 


FLO) BK Î 


Pe q+ AT V= KV)? (3-15) 


4 Po 


La figure (3-2) représente la densité f,(v). La fonction f(v) résulte d'une 
combinaison linéaire des courbes f, correspondant aux valeurs entières 
positives, négatives ou nulles de k avec en plus f,(v) = f,.(-v). 


IV. SCHÉMAS CONDUISANT A DES SPECTRES EN 1/9 DANS UN 
LARGE DOMAINE 


1°) Une remarque de M. SURDIN (14) 


Les schémas précédent donnent des spectres en 2 mais ne fournis- 
sent pas de région où la densité soit du type ee comme l'expérience le ré- 
clame. On peut cependant à partir des résultats précédents, obtenir des 
spectres en L dans un large domaine en utilisant une intéressante remar- 
que de M. SURDIN que nous reproduisons ici. Reprenons, par exemple, la 
relation (2-16) en posant 


+= (4-1) 


il vient 


FCV? œ? ie VE (4-2) 


Supposons alors, en employant le langage du schéma n° 1, que nous soyons 
en présence d'atomes non plus correspondant à une seule valeur de à mais 
à un ensemble continu de valeurs de à, tous ces atomes opèrant indépen- 
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damment. Soit g(x) dx le nombre d'atomes pour lesquels & est compris en- 
tre &x et à + dx. Le nouveau spectre aura une densité 


f'(v)= 2 [7 xg(x) da (4-3) 
U  G?+A4NR V? 
Supposons que l'on ait, 
PB: 
g[o=— 
Alors il vient 
fiv) tie sr tq (ie MR SE (4-4) 
Tv 2v O1 A; 
1+ AT V2 


Supposons de plus &, > 1. Etudions alors l'expression (4-4) dans les 3 cas 
particuliers suivants : 


a) 210 <<: 
Alors : fu(v):= 2 (4-5) 


la densité est constante. 


b) AX<2TV <a alors on a, tant que 2Fiv ne s'approche pas trop des ex- 
trémités de l'intervalle 


f'(v) s-È (4-6) 


2V 
on a bien ici une loi de variation en 1/24 


c) enfin si à,<@,< 2Tv, il vient 


FU) ER (4-7) 
2TAN2 
Une comparaison graphique des approximations (4-5) (4-6) (4-7) avec la 
courbe exacte (4-4) confirme d'ailleurs le bien fondé de ces affirmations. 
Il y a un domaine ou f'(v) varie à peu près en 1/v puis, pour les grandes 
valeurs dev , f'(v) tend à décroître comme 1/v:. 


2°) Les schémas de Mac-FARLANE (15) 


Cet auteur discute l'effet de scintillation des cathodes à oxydes. Dans 
une interprétation il considère un mécanisme du type suivant : l'émission 
d'une cathode à oxyde de baryum est sous la dépendance directe des atomes 
de baryum existant sur la surface émissive. En étudiant les fluctuations de 
la distribution de ces derniers qui résultent à la fois d'un effet de diffusion 
du baryum vers la surface et d'un effet de conduction ionique, il est con- 
duit à calculer une fonction de corrélation qui est très voisine de 


C0} 200 TPE EP (4-8) 
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où m est compris entre 0 et 0,5 et où k,{x] est la fonction de BESSEL mo- 
difiée de seconde espèce définie par 
(EN 9 ELLE [es GES 
2 sin Tv 
Ceci nous conduit à dire un mot des fonctions de corrélation du type (4-8) 


car, comme nous allons le voir, elles vont nous permettre de trouver, 
d'un seul coup, tous les spectres qui nous intéressent. 


Les fonctions K possèdent la propriété suivante : q et t étant positifs 
et m >-5 on a (Voir Wattson : Theory of Bessel functions). 


(4-9) 


2Ti vT 


TVA) TEE 


(gti) ka (gT)= VTT 2 (mehr f a (4-10) 


Le facteur 1/(q° + 41)" étant toujours positif, peut être interprété com- 
me une densité spectrale et, par conséquent, le premier nombre représen- 
te une fonction de corrélation. En faisant varier m, nous pouvons ainsi Ob- 
tenir un ensemble de paires de fonctions associées (l'une étant une densité 
et l'autre une fonction de corrélation). Le tableau ci-dessous précise quel- 
ques cas particuliers. Il donne, pour diverses valeurs de m, la fonction de 
corrélation c(T) et la densité spectrale f(v) correspondante. Pour m £< 0 
on a c(o)= +00. La puissance moyenne correspondante serait donc infinie, 
ce qui est physiquement inacceptable. 


PAIE a le spectre 
(q? + 47° v?) *? tend 


à devenir 
uniforme 


En fait, une limitation du spectre du côté des hautes fréquences (limitation 
qui existe dans tout appareil d'observation) fait disparaître cette difficulté. 


On peut donc dire que les fonctions correspondant à - D<m < 0 sont gfé 
aspects limites de fonctions physiquement acceptables. Lorsque m + - LÉ 
le spectre tend à devenir uniforme. Cela est naturel. En effet, pour 


CHR Le ik 
E 


m = + tandis que 
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est fini pour & petit, je [c(t)ldrest infini. Lorsque m + + la corréla- 
€ 


tion tend donc à devenir microscopique et le spectre évolue vers un spectre 
uniforme. Il faut, d'autre part, noter que 


Valti ke girl] =VT e°* (4-11) 


de sorte que le cas m = Lest identique à celui qui est décrit par les équa- 
tions (2-14) et (2-15). 


La figure (4-1) donne pour q = 1; T > 0 les courbes représentant les 
fonctions 


CN Ce TC (TA 


lorsque m prend l'une des valeurs 1, 1/2, 0. 


La figure (4-2) fournit, dans les mêmes conditions, les courbes repré- 
sentant f(V)/ f(0) pour v > 0. 


c(IS1) 


CONCLUSION 


Les schémas précédents n'abordent pas la question de l'explication 
physique précise de l'effet de scintillation; il nous a cependant paru utile 
de les rassembler afin de mettre en lumière divers mécanismes statisti- 
ques qui peuvent représenter, jusqu'à un certain point, les résultats expé- 
rimentaux. Nous pensons par ailleurs que ces modèles peuvent être utiles 
dans de nombreux problèmes présentés par l'étude de phénomènes de fluc- 
tuations dans des domaines très divers. 
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L'ESTIMATION DES PARAMÈTRES 
DES FONCTIONS D’ENGEL 


par 
Rémy FREIRE 


"In the case of qualities there was considerable evidence for a semi- 
logarithmic function which ..... leads to the equation 
(2 


pa; +b.1lo 
d'a La MS ken. 


LE 


In this case the error term might well be assumed to be normally and 
independently distributed but the estimation of the parameters is greatly 
complicated by non-linearity. Customary computation techniques are po- 
werless here and the only possibility seems to be an iterative solution 
using electronic equipment. " 


"The Econometrics of family budgets", H.S. HOUTHAKKER, Journal 
of the Royal Statistical Society, Series A, Part I, 1952, page 14. 


0. Le traitement statistique des fonctions d'Engel n'est pas un problème 
compliqué quand on accepte l'hypothèse simplifiée de la linéarité. Il est 
pourtant évident qu'une telle hypothèse ne peut être qu'une première appro- 
ximation, laquelle s'éloigne forcément de la réalité à mesure que le revenu 
augmente. Cette observation que l'on peut, d'ailleurs, accepter intuitive- 
ment sans peine, est largement confirmée par les données des études pra- 
tiques sur le sujet. 


En conséquence, pour avoir une image plus approchée de la réalité, on 
a essayé plusieurs expressions fonctionnelles non linéaires et, notamment, 
des polynomes et des fonctions logarithmiques. Bien que celles-ci appa- 
raissent dans les dernières études comme les plus adaptées au problème, 
il est possible de garder une certaine généralité en acceptant que les fonc- 
tions d'Engel, susceptibles de présenter un intérêt pratique, peuvent être 
“ésumées dans l'expression 


Yj= + B; RES TE, 


où x représente la dépense totale, y le prix spécifique pour le poste j du 
budget et, BA ...-xétant des paramètres caractéristiques, les x x 
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correspondant notamment à des coefficients d'équivalence entre les con- 
sommations des différentes catégories de personnes composant chaque 
grappe familiale. 


Cette amélioration théorique des fonctions d'Engel, bien que serrant 
de beaucoup plus près les données de la réalité économique, présente des 
difficultés d'ordre statistique dans l'estimation des paramètres et on est 
arrivé à envisager l'utilisation de nouvelles techniques de calcul par des 
machines électroniques comme étant les seules capables de porter un se- 
cours effectif à la résolution du problème. 


Dans ce travail il est précisément question de démontrer que les pro- 
cédés statistiques habituels,et notamment la méthode d'estimation par le 
maximum de vraisemblance, permettant le traitement du problème dans 
des conditions de calcul pratique parfaitement abordables. En plus, la mé- 
thode proposée conduit à une analyse détaillée de chaque poste spécifique 
du budget, ce qui permet de déceler et de traiter convenablement toutes 
les caractéristiques économiques qui composent le problème général. 


1. Si l'on admet que la fonction f(x,X,,X,, -..X,.) possède des dérivées 
secondes continues sur toute la région où la courbe d'Engel aura un sens 
économique et acceptant que les écarts entre les valeurs de la fonction en 
chaque point et celles fournies par l'observation se distribuent normale- 
ment et indépendamment, l'estimation par la méthode du maximum de 
vraisemblance correspond à chercher les valeurs des paramètres maximi- 
sant la somme : 


D [y-a-8 Flux.) |? 


dans laquelle la sommation se réfère à l'observation du même poste sur n 
budgets. C'est-à-dire, nous sommes conduits au système classique normal 
d'équations suivant où nous avons désigné par à, b, k,, k,, ... k,les in- 
connues représentant les estimations pour les paramètres «,8 ,4,,X, , 
... X_ respectivement : 


A1 io 22 +2y = 0 


—axf SOS tr VON 


of af af 
- 3b % 3 
tra El Er 


pit vnpé)hèt) Lo 


Suivant maintenant la technique classique, par dérivation successive 
de toutes les équations du système, en ordre à toutes les variables, chan- 
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geant tous ies signes et substituant à chaque vaieur de y sa valeur espérée, 
nous obtenons la matrice d'information : 


rte 


CLR of 
g Or. 

n DS D _ ba me 

d'A EP (a us 

DE DE boue De >= xs 

of LR NE TE of of 

URSS TER LI rs lof: 

ET > rot) DE el 


Par inversion, cette matrice fournit alors immédiatement la matrice 
des moments du deuxième ordre, c'est-à-dire, la matrice des covariances, 
que l'on peut écrire : 


(V) = [ Was Wab Wa Mars 
Wèb Wbb Wok Fe Wok my, 
W K4 Wu, Wie ti Mk, km/2 
7 Wok, Matmpe UT mime 


en symbolisant par W.., W, Wx,, ... W,.,les composantes de la matri- 
ce inverse de 


fn 5 né … rs 
Df DE E ee fit | 
of of orŸ of or 
Son an) (a fa 
of of [af af che 
= dk» Es EE ni Les (a) 


L'examen de cette matrice révèle un fait d'une importance centrale 
pour notre problème. En effet, ses composantes dépendent exclusivement 
de la fonction f et de ses dérivées, c'est-à-dire, la matrice (T) dépend seu- 
lement des valeurs des dépenses totales x et des estimations pour les para- 
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mètres X,, Xa, eee Xm» onc de x et de k;, k,, ...k,, étant indépendante de 
toute estimation pour les paramètres « et f. 
Soit maintenant a', b', k!, k!, ...k!,un ensemble d'estimations ineffi- 


= R— _ 

cientes pour tous les paramètres, obtenu par une méthode quelconque. 
Alors, si ces estimations ne sont pas exagérément grossières, en les in- 
troduisant dans les premiers membres des équations du système normal 
on obtiendra des sommes, lesquelles, bien que différentes de zéro, pré- 
senteront pourtant des valeurs relativement petites et que l'on peut dési- 
gner par A, , At, AK, AK, ee A,,, correspondant, dans cet ordre, à la 
première, deuxième, ... nième équation du système normal. 


Si l'on additionne alors aux estimations inefficientes a', b', k!, k’, 
Lt leurs incréments da: 0b 000 en -#10pue respectivement, on 0b- 
tient les vraies valeurs pour le maximum de vraisemblance à, 6, k,k, 
... K.. Ces valeurs n'étant pas connues, évidemment, il est nécessaire 
d'en obtenir une détermination approchée. C'est un des mérites de Fisher 
d'avoir montré qu'une telle détermination efficiente est donnée par les 


expressions : 


DAS NES AE IWES ALI TINET ER CREME REPARER Wars 


Wb JA. ETS Wbb NS ae Wu, AK, te Wbu, A4, A4 + WumAk 


©? 
Se 
Il 


ôk, = Wu A + W 


7: 40b + Wok du + W 


AR teatiN A! 


«ik Kikm4 — Km 


SRE Aa + Wok Db : Wu dk +W 


! 
lai ne AD test We x 


AS 
mb Km 


Maintenant si l'on considère les termes des équations du système nor- 
mal séparément, comme appartenant au cadre d'une matrice, on peut re- 
marquer que les deux premières colonnes de la matrice d'information mul- 
tipliées, respectivement, par (-a') et (-b'). En développant les expressions 
de À,, Ai, A,,A,, ... A, dans les formules pour 64,65 ,6, , ...8, , une 
partie se trouverà être égale au dernier terme de l'équation 'correspündan- 
te et, si d'accord avec la remarque sur l'égalité des colonnes du cadre du 
système normal et des colonnes correspondantes de la matrice d'informa- 
tion, on prend les S'6" comme solutions formelles d'un système d'équa- 
tions, leur expression devient : 


N'AURA _[, 2? 

Ga = 8 +WSEY+WeZ (JT) +U, na E(y a) 

ë, = -b'+WpEy+ My) +W Dire +...+W He 
5 pk, ok, bkm War 


)E. 
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ôk SE Era EN) ES + + W 
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à Comme les estimations efficientes pour & et B sont respectivement 
a" +0; et b' +6}, un ensemble complet d'estimations pour tous les para- 
mètres sera fourni par les expressions : 


a" — Was 2 Y + Wab2 (yf)+ wc ARTE M, 0 _… 
DAMES VE M2 (yf)+ Ma be PARLE Wbx, 2 b | 
et 
RDF MMS à 
SES ie VE Wu Z (YF) + WE Sa ue 

of 
mb] Gran 

Le TE 

"40 


L'inspection de ces expressions nous révèle le fait très important que 
les estimations efficientes sont fonctions exclusives des valeurs observées 
et des estimations préalables k', k;, ... k', pour les paramètres %,,x,, 
...Xn les estimations a' et b' pour & et B étant complètement absentes de 
ces formules. C'est,d'ailleurs, cette absence qui permet une application 
effective et pratique de la méthode, comme on va le vérifier par la suite. 
En effet, ce fût la généralisation d'une méthode basée sur une remarque 
semblable et proposée par le Prof. W.L. Stevens pour l'estimation des pa- 
ramètres de la régression exponentielle (Biometrics, Vol. 7, n° 3, 1951, 
"Asymptotic Regression") qui a permis d'envisager ce procédé pour l'esti- 
mation des paramètres des fonctions d'Engel. 


En conséquence, si par une certaine méthode on obtient un ensemble 
d'estimations non efficientes pour les paramètres X,,%X,, .:..Xm , On Cal- 
cule facilement un nouvel ensemble d'estimations efficientes pour tous les 
paramètres par le procédé particulièrement simple, lequel consiste à cal- 


culer les sommes 
of of 
2 Vue rs iTi) ARS EE 


et à multiplier successivement chacune d'elles par chaque composante W 
de la même colonne de la matrice[ V]et à additionner les produits. Utili- 
sant les deux premières colonnes de cette matrice on obtient des estima- 
tions efficientes pour les paramètres & et f8 ; les estimations efficientes 
pour les paramètres X,,X,, ... Xm résulteront, ensuite, par addition des 
sommes correspondant aux autres colonnes aux estimations inefficientes 
respectives, ceci après division par l'estimation trouvée pour le paramè- 
tre 8. Cela revient à dire que l'effort principal de calcul se concentre sur 
l'inversion de la matrice| T] qui fournira les composantes W de la matrice 
des covariances [V]|. | 

Si l'on part d'un ensemble d'estimations inefficientes très éloigné des 
valeurs du maximum de vraisemblance, il se peut que l'application du pro- 


cédé devienne nécessaire une deuxième fois. Pourtant, comme les estima- 
tions obtenues une première fois sont efficientes, cette deuxième approxi- 


mation sera sûrement suffisante. 
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Un indice de l'approximation des estimations obtenues est fourni par 
la première équation du système normal. En effet, si on est devant l'esti- 
mation du maximum de vraisemblance, la substitution de ces estimations 
dans le premier membre de cette équation doit fournir un résultat dont la 
différence à zéro doit être totalement attribuable aux erreurs d'arrondis- 
sement. Ceci, pourvu que tous les calculs soient corrects, ce qui veut dire 
que la même preuve peut servir aussi à la vérification de l'exactitude des 
calculs entrepris. 


Eñ résumant, on peut dire que les avantages de la nouvelle méthode 
proposée sont fondamentalement de deux sortes : du point de vue des cal- 
culs numériques, la question se réduit essentiellement à l'inversion d'une 
matrice, problème conventionnel pour la résolution duquel on possède plu- 
sieurs techniques parmi lesquelles le choix dépendra surtout des valeurs 
des composantes et de l'ordre de la matrice à invertir, du point de vue 
économique, comme le processus de convergence est indépendant des esti- 
mations pour «à et B, les résultats permettent d'isoler et de traiter conve- 
nablement le comportement de chacune des courbes d'Engel correspondant 
aux différents postes qui composent le budget et aussi l'analyse comparati- 
ve de leur évolution. 


2. Jusqu'ici on a envisagé une fonction d''Engel dans de très larges condi- 
tions de généralités. Pourtant, dans la pratique, les expressions analyti- 
ques concrètes utilisées sont presque exclusivement des polynomes et des 
fonctions logarithmiques et il apparaît que ces dernières sont peut-être les 
plus adaptées au problème. En particulier, si l'on reprend le type de fonc- 
tion utilisée par Houthakker dans son travail "The econometrics of family 
budgets", (Journal of the Royal Statistical Society, Series A, Part I, 1952, 
page 14) on peut écrire : 


Y= œx+P log Se 


. Utilisant les lettres a, b, k,, k,, ... k, pour symboliser les estima- 
tions des paramètres &, B,X,,%X,, +-. Xm respectivement et faisant 


1 
2piki 7% 
pour simplifier, on peut écrire les équations d'observation sous la forme 
commode 


y=a+b log (gx) 


Pour un premier ensemble d'estimations inefficientes des paramètres 
X49,X29, +++ Xm ON aura une première valeur g' deg et, pour le cas parti- 
culier en question, la matrice (T) prend la forme 


(T)=f on ZX log(g'x) -np,g' senps g! 
log (g'x) 3 log?(g'x) —p,g" Z log (9'x) .. .-p,9'Zlog(g%x) 
=np,g' -p,g'>log(g'x) np g'? + NP; Pm9'? 


En Pa SPAS 22 109 (dx) RP Ph EN Oe R0 pe dE 
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Invertissant cette matrice, on obtient une autre dont les composantes 
sont les éléments W de la matrice (V) des covariances, comme il fut dé- 
montré antérieurement. Alors en disposant en colonne, à gauche de cette 
dernière matrice, les éléments par lesquels il faut multiplier ses compo- 
santes, on peut disposer les calculs de la façon schématique suivante : 


Z y Was Web Wu, Sac Wars 

2 LD log (g'x)]| We Web Wok, Wu 

9 Zy Wa, Wok, Wk, DATE 
Pmg'Zy War, Wok, Wu 3 0ù Wet 


Multipliant successivement chaque élément de la colonne à gauche par 
chaque élément correspondant de la matrice et additionnant les produits on 
obtient, par l'entremise des deux premières colonnes de la matrice, des 
estimations efficientes a" et b" pour x et B respectivement. Des estima- 
tions efficientes k", k°, ... k® pour les paramètres %X,,X,, ... Xm résul- 
teront de la même façon, mais après division par l'estimation b" obtenue 
pour le paramètre. Œ 


Ensuite, la courbe d'Engel peut être graduée par la formule 
ÿe= a"+b" log(g"x) 
ou 
EE 
9 L 2 p: k': 


et les valeurs ainsi obtenues comparées aux données de l'observation. 


3. Comme déjà noté, si la substitution des estimations dans le premier 
membre de la première équation normale révèle une approximation insuf- 
fisante, on doit appliquer la méthode une nouvelle fois, mais en partant 
des estimations efficientes déjà obtenues. 


Une fois que les estimations obtenues soient suffisamment approchées, 
on peut passer au calcul de la variance résiduelle, que l'on obtient en divi- 
sant la somme des carrés des écarts 


Z(y.-y) 


par les n-m-2 degrés de liberté que comporte le problème. Par multiplica- 
tion par les composantes de la matrice { V] on aura, finalement toutes les 
variances et covariances des estimations trouvées. 


Ici se place, pourtant, une observation d'une certaine importance. En 
effet, la dernière matrice [V] effectivement calculée ne correspond pas 
exactement aux dernières estimations obtenues mais, précisément, aux 
avant-dernières et le calcul des moments du deuxième ordre de la façon 
indiquée ne conduit pas aux valeurs exactes de ces moments mais à des 
valeurs approchées. Mais, comme la loi de distribution pour les paramè- 
tres n'est pas connue, il sera nécessaire, en toute circonstance, de pren- 
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dre la loi du t' de "Student" comme distribution approchée ce qui nous con- 
duit forcément à une solution non exacte. C'est-à-dire, le travail de recal- 
culer la matrice [V] serait stérile, non seulement face à la nature du pro- 
blème des intervalles de confiance quine demande pas une solution exacte, 
mais aussi du fait qu'une solution précise ne peut pas être atteinte. 


Aussi, si l'on prend un paramètre quelconque, disons k;, les interval- 
les de confiance correspondant à sa détermination par l'entremise des dif- 
férents postes du budget, peuvent être comparés à un certain niveau de si- 
gnification - celui des 5 % , sauf motif justifiant un autre choix - et le com- 
portement relatif de chaque intervalle convenablement examiné.Théorique- 
ment, une région commune de recouvrement pour tous les paramètres de- 
vrait exister, dans le cas de compatibilité de toutes les déterminations 
mais, la nature même des différents postes du budget doit conduire, prati- 
quement, à une situation où une telle éventualité ne se présentera pas. Ceci 
veut dire que l'interprétation des écarts que certaines des courbes présen- 
teront face à la région de recouvrement et aussi le choix de la meilleure 
valeur de compromis pour les différentes estimations de chaque paramètre 
sera donc un problème plus économique que statistique. 


Finalement, si l'on prétend obtenir un ensemble d'estimations pour 
tous les paramètres et correspondant à l'ensemble de toutes les courbes 
d'Engel, on peut, après avoir fixé les estimations trouvées comme les 
plus convenables pour les paramètres X,, X,, +... Xm, reprendre une der- 
nière fois le problème de l'estimation des paramètres « et 8 par une sim- 
ple technique d'ajustement linéaire. 


ÉTUDE STATISTIQUE DE CONTOURS CRANIENS 
CONSIDERES COMME COURBES ALÉATOIRES 


par 


Marie-Jeanne LAURENT-DUHAMEL 


INT RODUC TION 


Dans son étude générale sur les éléments aléatoires de nature quel- 
conque, M. le Professeur FRECHET a envisagé l'étude de courbes aléa- 
toires planes fermées. Sous sa direction ont été étudiées, en particulier, 
deux séries de 100 courbes obtenues en prenant à l'aide d'un conformateur 
de chapelier, d'une part les contours de têtes de 100 hommes adultes au 
niveau habituel des chapeliers, d'autre part les contours de 100 crânes 
humains, pris dans une collection du Musée de l'Homme et provenant 
d'un ossuaire des Iles Canaries. Ces contours ont été relevés dans un 
plan horizontal passant au niveau supérieur de la soudure des parié- 
taux et des temporaux, (la verticale étant définie par la droite quijoint 
les incisives supérieures au point le plus saillant de l'os frontal ). 
Bien que ces dernières courbes ne présentent peut-être pas une pré- 
cision de détermination aussi rigoureuse qu'on pourrait le souhai- 
ter du point de vue anatomique , elle présentent sur les précéden - 
tes une certaine supériorité d'homogénéité en raison du choix du 
metériel de base et de la manipulation continue du conformateur 


—{C) : contour donne par 
l'appareil, 
(09: contour ree/ 
H reconstitué 
+ Q : centre de | appareil 
! pôle de là transtor 
mation 
AA: dxe de l'apparell 
pris comme axe polaire. 


par le même manipulateur pour le relevé de ces contours. L'appareil, dé- 
nommé conformateur, ne donne pas directement les contours de têtes, 
mais des contours réduits et déformés obtenus en retranchant de chaque 
rayon passant par un point fixe de l'appareil, repéré à l'intérieur des cour- 
bes, un segment de longueur constante, égale à celle des branches mobiles 
de l'appareil, soit par exemple pour celui qui a servi au relevé des con- 
tours crâniens, 54 mm (fig. 1). On obtient le contour réel par la transfor- 
mation inverse connue en géométrie sous le nom de transformation con- 
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choidale. On notera que cette transformation peut changer assez fortement 
l'allure de la courbe en accentuant ou en créant des irrégularités telles 
que : changement de concavité, points de rebroussement,... Singularités 
bien connues des conchoïdes classiques (limaçon de PASCAL, cardioide. :) 
de la plus régulière des courbes : le cercle. 


Il. ÉTUDE MATHÉMATIQUE DES CONTOURS 


Pour caractériser analytiquement ces courbes, il a paru simple du 
point de vue mathématique, de former pour chacune son équation en coor- 
données polaires P =f (w). On a pris pour pôle le centre de l'appareil, 
point défini plus haut comme pôle de la transformation conchoïdale, immé- 
diatement repérable sur les courbes, et pour axe polaire, l'axe de symé- 
trie de l'appareil au repos, lui aussi directement repérable sur les cour- 
bes; l'avantage immédiat de ce système de référence résidant dans le pas- 
sage très simple de l'équation d'un contour transformé donné par l'appareil 
à l'équation du contour vrai. 


Pour obtenir une fonction analytique représentant chaque courbe, on a 
calculé, en prenant 24 points sur chaque contour et les rayons vecteurs 
correspondants, les coefficients des développements en séries trigonomé- 
triques des fonctions périodiques P (w) de période 27 . Ces développements 
ont été arrêtés au terme d'ordre 12, les coefficients d'ordres supérieurs 
étant généralement négligeables devant les autres. On a ainsi obtenu pour 
chaque contour 24 coefficients; les coordonnées d'un point (P;:w) du contour 
jG=1, 2, 3, ... 100) étant liées par la relation : 


\1 
P;= 24 œ;, cosi w+2 b;; sin; w 
nn nb ec 


En raison du choix du pôle, seul le terme en a, varie (et d'une quanti- 
té constante, 54 mm) quand on passe de la courbe brute au contour réel par 
la transformation conchoïdale. 


Les courbes brutes (transformées des contours réels de têtes et de 
crânes) ont été microfilmées, ainsi que les 100 valeurs des 24 coefficients 
des développements trigonométriques, pour les deux séries de courbes. La 
reproduction de ces courbes et des tableaux des coefficients étant trop 
étendue pour figurer dans cet exposé, les microfilms ont été déposés au 
Laboratoire de Calcul de l'Institut HENRI POINCARE où ils peuvent être 
consultés ou donner lieu à des reproductions photographiques. D'autres 
courbes aléatoires, ét:diées par M. le Professeur FRECHET, obtenues 
en jetant 100 fois un fil fermé sur une surface plane, ont été également 
microfilmées ainsi que les courbes donnant les variations des coordonnées 
cartésiennes d'un point de chacune de ces courbes. Les développements en 
séries trigonométriques de ces dernières ont été entrepris. Les micro- 
films ont été joints aux précédents. 


II. ÉTUDE STATISTIQUE DES COEFFICIENTS DE DÉVELOPPE- 
MENTS 


Chaque coefficient se présente comme une variable aléatoire dont nous 
possédons cent valeurs. Les lois de distribution de chacun de ces coeffi- 
cients peuvent-elles être considérées comme normales? Pour répondre à 
cette question nous avons utilisé le test du X° et le test de normalité de 
GEARY et PEARSON. Les résultats figurent aux tableaux I et Il. 
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Les conclusions que l'on peut en tirer concordent assez bienet sem- 
blent permettre de répondre par l'affirmative sauf peut-être pour le terme 
b, pour lequel, dans le cas des contours de têtes, le test du X° donne un 
niveau de probabilité extrêmement faible, alors que le test de GEARY don- 
ne des niveaux de probabilités acceptables; dans le cas des contours crâ- 
niens, les niveaux de probabilités pour b donnés par les deux tests sont 
faibles. A part cette exception qui porte sur le dernier coefficient des dé- 
veloppements, assez petit en général pour que l'on puisse négliger le ter- 
me correspondant, on peut admettre que les coefficients des développe- 
ments en série trigonométrique suivent des lois de répartition normales. 


Dévendance des coefficients 


Parmi les 276 corrélations existant entre les 24 coefficients pris deux 
à deux, nous en avons étudié une dizaine pour chacune des deux séries de 
courbes. Pour cette étude, on a calculé les coefficients de corrélation, les 
rapports de corrélation, les indices diagonaux de M. FRECHET, les indi- 
ces de connection de GINI, les indices de JORDAN et les indices de 
M. GEIRINGER (Tableau Ill). 

Les résultats obtenus ne permettent pas de conclure à une dépendance 
très forte entre les coefficients des développements. On peut remarquer 
que les quatre derniers indices varient relativement peu, et restent du 
même ordre de grandeur. 


IV. ÉTUDE DES CONTOURS CRANIENS CONSIDÉRÉS COMME 
COURBES ALÉATOIRES 


Nous avons essayé de définir des contours typiques analogues aux va- 
leurs typiques classiques des variables aléatoires. Soient : un contour 
moyen, un contour médian et des contours quartils. 


Définition d'un contour moyen 


Les coordonnées polaires étant celles précédemment définies pour 
chaque contour, nous conviendrons d'appeler "contour moyen", un contour 
tel qu'à partir d'un axe et d'un pôle donnés, on a pour chaque valeur de 
l'angle polaire, le rayon vecteur correspondant de ce contour moyen : 

100 


p (w)= 2 P; (uw) 


P;(w) étant le rayon vecteur correspondant à l'angle w du contour j. 


Cette définition du contour moyen ne s'impose nullement et pourrait 
être différente (ne serait-ce, par exemple, qu'en échangeant le pôle des 
coordonnées) mais elle semble raisonnable (voir V, paragraphes 1 et 2). 

On obtient immédiatement les coefficients du développement en série 
trigonométrique de p comme moyennes des coefficients de même rang de 
tous les développements : 


a; et b;; étant les coefficients de rang i dans le développement du con- 
tour j. 
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Connaissant les coefficients du développement dep , on a pu reconsti- 
tuer le contour moyen et obtenir le tracé, à partir de 24 points, soit de la 
transformée, soit du contour en vraie grandeur. 


Les contours médians et quartils sont définis de manière analogue; 
(w) par exemple, étant la valeur médiane des différentes valeurs de p; 
pour un angle polaire (w); les coefficients de rang i du développement de 
sont les valeurs médianes des coefficients de rang i des développements 
des p;. 


De même pour les quartils. On a ainsi reconstitué et tracé quatre con- 
tours typiques. Le premier mode de comparaison qui semble s'imposer en 
raison de la manière même dont ils ont été obtenus, a été de superposer 
ces quatre contours en faisant coïncider les axes polaires et les pôles (fig. 
2). On constate que : 


1°)- Les Contours moyen et médian diffèrent très peu l'un de l'autre, 
ainsi que pouvait le laisser prévoir le fait que les coefficients a et b ayant 
des distributions sensiblement normales, leurs moyennes et leurs média- 
nes diffèrent très peu. | 


contour moyen 
a 1er quartil 
——. jène quartil 
contour Mmedian QU! 
coïncide , sauf entre contour moyen 
A el B,avec le 19 quartil. 
contour moyen. ——— 3tre quertil. 


Fig. 2-Comparaison des contours typiques Fig.3-Comparaison des contours typiques 
par coïncidence des pôles et de l'axepolaire. par coïncidence de l'axe de symétrie et de 
la binormale aux contours. 


2°)-Les quatre contours typiques semblent coïncider assez bien du 
côté de l'occipital et du temporal droit; les différences étant plus sensibles 
du côté du temporal gauche et surtout du côté du frontal. Ceci s'explique 
par la simple observation du procédé de relevé des contours et du choix 
des axes de référence. En effet, la différence inconsciente des rôles joués 
par les mains d'un manipulateur droitier dans le relevé d'un contour, dé- 
truit la symétrie de l'appareil au repos. Si les axes de référence sont liés 


à l'appareil, ce qui est le cas, et non au crâne, on doit retrouver cette dis- 
symétrie dans la comparaison des contours quand on superpose ces axes 
de référence. 


Une étude ultérieure d'ajustement de ces contours par des courbes dé- 
finies géométriquement nous a conduits à rapporter les contours à des axes 
définis, d'une part, par l'axe de symétrie, qui existe toujours très sensi- 
blement dans chaque contour et, d'autre part, par la perpendiculaire à cet 
axe menée par les points où les tangentes au contour sont parallèles à l'axe 
de symétrie. La comparaison des contours typiques rapportés à ces axes 
est donnée par la figure (3). 
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V. ÉTUDE DE LA VARIATION DU CONTOUR MOYEN AVEC LE 
CHOIX DES AXES DE RÉFÉRENCES 


1. Un changement d'axe polaire qui conserve le pôle ne change pas 
le contour moyen. 


Pour un rayon vecteur donné, le centre de gravité des points des diffé- 
rents contours situés sur ce rayon vecteur est indépendant de l'angle polai- 
re. Or, par sa définition même, le point du contour moyen situé sur ce ra- 
yon vecteur coïncide avec le centre de gravité. 


2. Variation du contour moyen avec la position du pôle. 

On démontre que si M est un point du contour moyen correspondant au 
pôle C et si C' est un pôle voisin de O (c'est-à-dire si OO" est petit devant 
OM) M est, en première approximation, un point du contour moyen corres- 
pondant au pôle O' (c'est-à-dire que la distance MM! est un infiniment petit 
du deuxième ordre par rapport à OO"). On peut donc considérer le crâne 
moyen précédemment défini comme assez stable quand les coordonnées 
varient peu. 


VI. COMPARAISON D'UN CONTOUR QUELCONQUE AVEC LE 
CONTOUR MOYEN 


Pour une valeur donnée w de l'angle polaire,p; et P étant les rayons 
vecteurs du contour j et du contour moyen, la différence p; -P donne pour 
chaque valeur de w et pour chaque contour la distance polaire entre un point 
de ce contour et le point correspondant du contour moyen. L'intégrale 


v (P; - P) du donnerait la somme algébrique de ces distances le long des 


deux courbes mais pourrait être nulle pour deux contours très différents. 
Nous avons calculé pour chaque contour (j) : 


= J (P-P) d'u 


quine s'annule que pour P; = P quel que soit w et qui est égal à l? si P;-P =1, 
où 1 = constante indépendante de w. 


REMARQUE.- X° est invariant dans la transformation conchoïdale ayant 
pour centre le pôle des coordonnées, les calculs peuvent donc être effectués 
sur les courbes réduites et les résultats appliqués aux courbes vraies. 


Pour calculer Ke: nous remarquons en développant son expression que : 


; a 2H 
= Ls+ts- À) PPdu 


2 
Tr 


en appelant S; et S les surfaces comprises respectivement dans le contour j 
et dans le contour moyen. On obtient une valeur approchée de l'intégrale 


{ 2T ps EN = £ 
Je FA ae ET 22 2 a; +22 (85 8 +b; b) 


PE 
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à l'aide des coefficients du développement. On obtient une valeur approchée 
de S; soit en planimétrant la courbe j, soit en la calculant à l'aide des coef- 
ficients du développement (tableau IV) les différences relatives dépassant 
rarement 1/100 ce qui permet en particulier de justifier l'arrêt à l'ordre 12 
des développements en séries. On prendra pour S; l'une ou l'autre des deux 
valeurs ou mieux leur moyenne arithmétique. 


La valeur moyenne des X° s'exprime simplement en fonction de la sur- 
face du contour moyen S et de la moyenne des surfaces soit S. On trouve en 
effet en sommant les X? que la somme des termes en J{"P;Pdw , s'expri- 
me simplement en fonction de la surface S et l'on obtient 


va lys 
= Xi=(S-S) 


X° ne dépend donc des axes de référence que dans la mesure où S en dépend, 
et nous avons vu que si l'origine varie peu le contour moyen tel qu'il a été 
défini est assez stable. 


Si nous considérons les valeurs de X° comme servent à caractériser 
les écarts entre un contour j et le contour moyen, nous pouvons étudier la 
répartition des X, = VX2. Nous avons, à cet effet, calculé la moyenne et l'é- 
cart type des X; (tableau V). L'application du test du X° directement à la ré- 
partition des X ; donne un résultat qui permettrait de conclure à lanormalité, 
mais en raison du fait qu'on ne peut à priori affecter les X; d'un signe posi- 
tif ou négatif, il a semblé préférable de tenter la normalité non des X; mais 
celle des log X;. On trouve, en les groupant en dix classes de fréquences 
constantes, une valeur de X° = 3,52 qui, pour 9 degrés de liberté, corres- 
pond à un niveau de probabilité de l'ordre de 0, 93. 
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Tableau |  Contours de têtes 
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Tableau lil 


Indice 
de 


Note : le test de Fisher appliqué à la linéarité de la régression pour chaque 
corrélation étudiée donne un niveau de probabilité supérieur à 5 %. 
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Tableau IV 
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Tableau V 
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SUR L'AMORTISSEMENT DES FLUCTUATIONS 
ECONOMIQUES 


(Examen critique de la théorie de Tinbergen)® 


par 


G. MALÉCOT 


L'économiste hollandais Tingerben s'est efforcé, dans ses nombreux 
travaux (1), de mener l'analyse des fluctuations de l'activité économique 
d'une façon semi-empirique, Semi théorique. 

L'échec des tentatives purement empiriques comme celle de Harvard- 
incapable en 1929 de déceler la crise naissante - a été flagrant. 


U Quant aux schémas purement théoriques de l'économique rationnelle, 
ils reposent souvent sur des simplifications excessives qui rendent toute 


prévision correcte impossible. 


L'originalité de Tinbergen est d'avoir essayé de faire, dans l'évolu- 
tion économique, la part des circonstances imprévisibles ou inanalysables, 
avec les perturbations "aléatoires" qu'elles engendrent, et la part des en- 
chaînements fonctionnels et déterministes qui traduiraient les "lois écono- 
miques" d'une époque donnée, enchaînements que l'on demanderait à l'ana- 
lyse statistique de dégager avec précision en décomposant la variation de 
chaque grandeur économique à expliquer en une variation causale due à 
l'influence des autres grandeurs considérées, et un résidu aléatoire tra- 
duisant l'effet sur l'économie des chocs fortuits ou des causes négligées. 
La forme dela décomposition est certes inspirée par les conceptions à prio- 
ri de la science économique qui fixent le nombre et la nature des relations 
causales (relations de structure), mais la détermination précise des va- 
leurs numériques à une époque donnée des coefficients qui figurent dans 
ces relations se fait par des méthodes statistiques telles que la méthode de 
corrélation multiple. Il y a certes une part d'arbitraire, Tinbergen le re- 
connaît, dans cette séparation statistique de la partie causale et de la par- 
tie aléatoire, en particulier quand les chocs aléatoires se répètent. L'ana- 
lyse statistique n'aura quelque valeur que si les chocs sont assez espacés 
ou ont des effets d'assez courte durée pour que les réactions internes de 
l'économie, à la suite de chaque choc, puissent être observées et mesu- 
rées de façon adéquate. 

Les délicats problèmes de la mesure et de l'estimation des variables 
économiques ont fait l'objet de nombreuses recherches de Frisch et 
d'Haavelmo, et de la "Cowles Commission for Research in Economics " 


(Lawrence Klein) (2). 


(a) Conférence faite à la Société Lyonnaise d'Econométrie le 25 mars 1953. 
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Je me contenterai, sans entrer dans des détails statistiques, d'exa- 
miner le modèle mathématique que Tinbergen a construit pour représen- 
ter le fonctionnement de l'économie des Etats-Unis entre 1919 et 1932 
(modèle que Lawrence Klein a remanié en tenant compte des données éco- 
nomiques jusqu'en 1941, mais dont la valeur de prévision pour l'après- 
guerre est naturellement reconnue extrêmement douteuse, même si on le 
complète, comme l'a fait Christ, .en tenant compte des données économi- 
ques de 1946 et 1947). 


Le point de départ dans la construction du modèle est la distinction 
entre, d'une part, les variables économiques dont il s'agit d'"expliquer" 
les fluctuations (production, prix, Salaires, etc...), variables que Tin- 
bergen appelle "endogènes", et les données non économiques (politiques, 
institutionnelles, psychologiques) que Tinbergen appelle extérieures au 
système ou "exogènes" (il ya là une distinction analogue à celle des forces 
intérieures et des forces extérieures de la mécanique). L'état et les va- 
riations des variables exogènes seront, comme le sont souvent les forces 
extérieures de la mécanique, regardés comme donnés et cela en général 
indépendamment des valeurs des variables endogènes, ces dernières étant 
définies par des "équations de structure" qui expriment comment les va- 
riables endogènes sont déterminées à partir des variables exogènes. 


Si l'on désigne les variables endogènes par x,, ...... , Xp et les va- 
riables exogènes par y, , ...... » Ya, 1eS équations de structure sont de la 
forme : 

| É (2 perse Xp» Y nee... Ya)= 0 
(1) f, COR RE rXprYire ess »Ya)=0 


et elles devront être en nombre égal à celui (p) des variables endogènes si 
l'on veut pouvoir calculer toutes celles-ci. La forme de ces équations sera 
donnée par la théorie économique, mais les valeurs numériques de cer- 
tains des coefficients (des élasticités en particulier) qui y figurent ne pour- 
ront être données que par l'observation statistique (avec toutes les limita- 
tions que celle-ci comporte, comme nous le soulignerons plus tard). 


Si les équations (1) ne reliaient que les valeurs, à un même instant t 
de toutes les variables, elles pourraient être résolues algébriquement en 
fournissant x,(t) ...... x(t) en fonction des variables exogènes PE) na 
+.., ya(t) Supposées connues. Plaçons-nous dans le cas particulier où les 
variables exogènes (institutions politiques, développement techni ue, be- 
soins des consommateurs et psychologie des producteurs, etc...) reste- 
raient constantes au cours du temps (ni mouvement séculaire, ni mouve- 
ment saisonnier), les variables endogènes seraient, elles aussi, constan- 
tes et définiraient l'équilibre économique du système. Mais on sait bien 
qu'une telle constance absolue est impossible à rencontrer dans la prati- 
que : les phénomènes économiques sont influencés, à côté des grands fac- 
teurs exogènes y;, par de nombreux facteurs secondaires que l'on ne peut 
faire entrer dans les équations (1) parce qu'ils sont trop nombreux, parce 
que l'influence de chacun d'eux est trop faible pour être pratiquement dé- 
terminable, et parce que leur sens de variation au cours du temps est im- 
prévisible; c'est à l'ensemble de ces facteurs que l'on donne le nom de 
“perturbations aléatoires", pour souligner qu'ils sont regardés, à chaque 
époque t, comme dépendant du hasard (exemple : fluctuation dans le temps 
de la demande par les consommateurs de tel ou tel produit) 
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Il est possible de corriger les équations (1) pour tenir compte des per- 
turbations aléatoires. Il suffit de les écrire : 


EE ne A nt, 


Fu » en D'ÉRNERE ÿ= up (t) 


Les u;(t) étant des variables aléatoires ou "chocs" qui traduisent les 
perturbations exercées par les causes aléatoires sur les p équations de 
structure. Les hypothèses que l'on fait en économétrie sur les u;(t) sont en 
général assez arbitraires. On suppose souvent (mais les vérifications sta- 
tistiques à l'appui sont peu nombreuses : voir toutefois l'ouvrage de 
Klein (2)que ni 9 et u;(t) sont stochastiquement indépendantes, ainsi que 
ui(t)etu.(t') (t'£t), et même que les u;(t) obéissent à des lois de Causs 
que l'on peut toujours supposer de moyenne zéro, car cela revient à faire 
rentrer dans les f; une constante convenable); on sait que l'hypothèse gaus- 
sienne est admissible dans la mesure où les influences aléatoires qui agis- 
sent sont additives et nombreuses, la part de chacune d'elles étant faible 
par rapport à l'ensemble. 


L'introduction dans les équations (2) des u:(t) entraïînerait évidemment, 
même si les variables exogènes y; étaient constantes, des fluctuations des 
x; au cours du temps, autour de l'état d'équilibre. 


L'étude de ces fluctuations serait du ressort de la théorie des "fonc- 
tions aléatoires"; mais ce n'est pas sur elle que s'était portée jusqu'ici 
l'attention des économétriciens (et cela tient évidemment à ce qu'on ne 
sait pas grand chose sur les u;(t)). 


Tinbergen, en particulier, néglige d'écrire les termes aléatoires 
ui(t), et raisonne donc sur les équations de la forme (1), c'est-à-dire sur 
des équations "déterministes" (la statistique n'intervient, dans son mode 
de présentation, que pour le calcul numérique des coefficients qui figurent 
dans ces équations). Il estime que les fluctuations économiques que l'ob- 
servation constate - même dans une période où les variables exogènes sont 
sensiblement constantes - résultent avant tout, non du hasard, mais du fait 
que certaines des variables endogènes à l'époque t sont influencées par les 
valeurs d'autres variables endogènes à des époques antérieures à t (par 
exemple, dans les budgets rion ouvriers, les dépenses actuelles dépendent 
plus fortement du revenu touché il y a quelques mois que du revenu actuel). 


Les équations (1), quand on les écrira explicitement, feront alors in- 
tervenir, outre les valeurs "instantanées" x;(t), les valeurs des x; à des 
époques antérieures à t, soient x; (t-œx), x; (t-fB), etc..., les constantes 
œ, B, etc... désignant les "décalages dans le temps" ou "time-lags", 
quantités qui ne peuvent d'ailleurs en général qu'être déterminées approxi- 
mativement par l'analyse statistique. 


Les équations (1) seront alors, non pas des équations algébriques (au 
sens élémentaire du mot), mais des équations aux différences finies que 
l'on pourrait d'ailleurs sans inconvénient - avec toute approximation fixée 
d'avance - remplacer par des équations différentielles (le remplacement 
des différences finies par des dérivées et inversement est d'usage courant 


en mathématiques appliquées). 
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Supposons que les équations (1) s'écrivent : 


PEL MORE ou) GE A IT Etc 


(3) ......e So ee 5 es seine sise ss 000 .e eee eee s els;r.ejelele eee ee e 000 


ee 
Se ee ee 0 0 0 s ou c'e 0e ee eo 0e e © eo vle e/e 0 eue) else 0 [e;0 ee 5070009 


(avec des seconds membres u;(t), on aurait des "équations aux différences 
finies stochastiques" ou des "séries temporelles autorégressives ; mais 
Tinbergen, nous l'avons déjà dit, fait abstraction des ui (t). 


Ces équations (équations "dynamiques" au sens de Ragnar Frisch) dif- 
fèrent des équations "statiques" telles que 


ÊLE RODPOER CE RESCHEE ROLE EL À ]—=0 


en ce que, même si l'on suppose que les variables exogènes y;(t) restent 
constantes au cours du temps, les x;(t), définis par (3), ne prennent pas 
instantanément des valeurs d'équilibre; si leurs valeurs initiales à l'ins- 
tant t = 0 ne sont pas des valeurs d'équilibre, les x;(t) peuvent, ou bien 
s'éloigner de plus en plus des valeurs d'équilibre (le système est alors dit 
instable) ou, au contraire, s'en rapprocher indéfiniment soit par un mouve- 
ment de sens constant, soit par un mouvement oscillatoire (le système est 
alors dit stable). (D'autres circonstances sont mathématiquement possibles 
mais nous les laisserons de côté comme étant sans intérêt pratique). 


Le problème que s'est posé Tinbergen est le suivant : déterminer sur 
les équations de structure, construites à partir des données économiques 
des U.S.A. entre 1919 et 1932, si elles auraient défini, en l'absence de 
toute variation des variables exogènes (pas de variation des y;) et de toute 
perturbation aléatoire (u; tous nuls), un système stable, et étudier dans 
cette éventualité les petits mouvements au voisinage de leur position d'é 
quilibre, en particulier du point de vue de leur facteur d'amortissement. 


Si leur amortissement est rapide, on peut penser - quoique cela reste en- 
core à préciser mathématiquement -[en dehors du cas linéaire étudié par 
Mann et Wald (Econometrica 1943.)](Voir aussi Tinbergen page 33) - 
qu'en général des variations exogènes et des perturbations aléatoires suf- 
fisamment faibles ou suffisamment espacées seraient bien ‘"encaissées" 
par le système qui parviendrait assez vite à en amortir les effets (en ex- 


cluant le cas évident d'oscillations forcées provoquées par des variations 
exogènes de période égale à une période propre des mouvements endogènes). 


Or, il est bien connu que, pour étudier la stabilité d'une position 
d'équilibre ainsi que les petits mouvements dans son voisinage, on rempla- 
ce les équations données par des équations approchées linéaires que doi- 
vent vérifier sensiblement les écarts des variables par rapport à leurs 
valeurs d'équilibre, pourvu que ces écarts soient süffisamment petits. 
Supposons donc déterminées les valeurs d'équilibre x: des x;(t) correspon- 
dant à des variables exogènes y; constantes données (en fait Tinbergen se 
contente d'une détermination statistique, en prenant pour valeurs d'équili- 
bre des variables endogènes x;(t) les moyennes arithmétiques de leurs va- 
leurs numériques pendant les 14 années de 1919 à 1932; si la compensation, 
sur ces 14 années, des petites perturbations déterministes ou aléatoires 
est admissible, il apparaît tout à fait illégitime de négliger l'effet sur la 
moyenne de la grande perturbation de 1929 et cela quelles que soient les 
causes, tant discutées par les économistes, de cette perturbation) 
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Faisons abstraction des difficultés pratiques et supposons connues les 
valeurs d'équilibre x;. On peut alors poser dans les formules (3) 
Ait) = xx Ules x, désignant les écarts à la date t par rapport aux va- 
leurs d'équilibre x;. 


En remplaçant ensuite les premiers membres des équations (3) par 
leurs développements limités du premier ordre par rapport aux x;(t) et 
tenant compte de ce que (par définition des valeurs d'équilibre) : 


PART tn DR TR 2e ]=0 


000000000000 


RER 


on obtient le système approché linéarisé (valable pour les x; petits, c'est- 
à-dire pour les petits écarts par rapport aux valeurs d'équilibre) : 


À Xe Àoat-ot À3 X,t-ptess..e = 0 


On a là un système linéaire d'équations aux différences finies (trans- 
formable, si l'on préfère, en un système différentiel). 


Ce système est linéaire à coefficients constants, et homogène; sa ré- 
solution est facile : la détermination des solutions, réelles ou imaginaires, 
de son équation caractéristique fournit comme intégrale générale une com- 
binaison linéaire de fonctions exponentielles, réelles ou imaginaires, du 
temps. 


La condition de stabilité est que toutes ces exponentielles tendent vers 
zéro quand t tend vers +c, c'est-à-dire qu'elles aient des modules < 1. (Si 
les exponentielles sont imaginaires, le mouvement correspondant admet 
des périodes ou des pseudo-périodes dues aux facteurs sinusoiïdaux). 


Pour les U.S.A. de 1919 à 1932, les équations linéarisées utilisées 
par Tinbergen sont les suivantes (il importe de ne pas oublier que ces 
équations ne définissent que les écarts des variables endogènes par rap- 
port à leurs valeurs d'équilibre, mais non les valeurs d'équilibre de ces 
variables, ni les "variables exogènes" qui sont - abusivement, peut-être, 
en ce qui concerne certaines des équations - regardées comme constantes) : 

Désignons par at) = à + a, le nombre d'heures de travail réalisées par 
unité de temps dans le pays, à l'époque t, nombre qui se décompose en 
b + b, heures relatives à la production de biens d'investissement nouveaux 
et © + c, heures relatives à la production de biens de consommation et au 
remplacement des biens usagés. 


Si 1 + 1, est le salaire horaire d'un ouvrier et L + L, le montant total 
des salaires versée, par unité de temps, on a évidemment : 


a+a,=b+b,+C+c,, d'où (puisque à = b+c): 


(6) 2 as=b,;+c 


(a) Nous conservons la numérotation et les notations de Tinbergen. L'unité de temps est 
prise égale à 4 mois. 
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Tinbergen a choisi ses unités (1) de façon que a=1 et 1 = 1; d'où : 
(1) L,=1l,#a, 


Si l'on désigne par ü + u, la quantité de biens de consommation vendus par 
unité de temps, à l'époque t, par U + U, leur valeur totale et par P + p. 
leur prix unitaire (qui coïncide avec l'indice du coût de la vie), on a : 


U+ U:=(P + pe) (ü + ui) 


d'où par linéarisation (les unités étant choisies de façon que p = 2'et 
u = 0,9) : 


(3) U,= 0,9p,+2u, 


Tinbergen admet ensuite que les ouvriers consacrent à des biens de 
consommation immédiate la totalité de leur salaire, alors que les sujets 
non ouvriers, qui reçoivent par unité de temps un revenu E + E;, en épar- 
gnent une partie E" +E! pour ne dépenser à des biens de consommation 
que E'+E!, avec d'ailleurs un décalage de huit mois environ. D'où l'iden- 
tité : 


(2) Etembipes 
et l'équation de consommation : 
(5) DElrr > 


Il est à remarquer que l'équation (5) que nous venons d'écrire est la 
première qui ne soit pas vérifiée rigoureusement, c'est-à-dire ne soit pas 
une "identité", une "relation comptable". C'est donc la première dans la- 
quelle nous aurions, dans une théorie plus raffinée que celle de Tinbergen, 
à introduire une perturbation aléatoire traduisant l'étalement au hasard la 
la consommation au cours du temps. 


L'équation (5) est aussi la première qui fasse intervenir un décalage 
dans le temps (on concevrait difficilement d'ailleurs, dans des relations 
comptables, de véritables décalages dans le temps). 

Tinbergen admet encore que le nombre d'heures de travail consacré 
(par unité de temps) à la production de biens de consommation est propor- 
tionnel à la quantité produite, c'est-à-dire que : 


C + ct=k(u + un) 


comme il a pris u = 0,9 et qu'il a calculé € = 0,8 (c'est-à-dire 80 % du 
nombre d'heures total, puisque a = 1), onak=0,9, d'où : 


(7) ct = 0,9 u, 
En outre, les écarts p, des prix à la moyenne sont regardés - sur la 


base, d'ailleurs médiocre, des statistiques utilisées - comme proportion- 
nels à la variation des quantités vendues, avec un certain décalage : 


(£) Pt = Te (rt >0) 
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On aurait fortement tendance à discuter cette équation. Tinbergen, lui- 
même, reconnaît qu'elle est fort approximative; quand le prix augmente, 
l'offre augmente mais la demande devrait diminuer (1) et le prix se réta- 
blir (c'est la théorie classique du prix d'équilibre); il faut donc admettre 
que les courbes de demande se déplacent par translation, alors que les 
courbes d'offre restent relativement fixes, ce que Tinbergen justifie dans 
le cas où il y a tendance monopoliste. Tinbergen remplace d'ailleurs la 
notion de "courbe d'offre" par celle de "courbe de fixation du prix" (Eco- 
nometrics, p. 32) et il admet que les prix sont fixés en fonction croissante 
(peut-être en vertu de la loi des rendements décroissants) du montant des 
PC CE passées, ce qui explique le décalage introduit dans l'équation 


L'équation suivante détermine le salaire horaire des ouvriers, 1 + 1, 
en fonction du coût de la vie p + pet de "l'emploi", c'est-à-dire du nombre 
total d'heures de travail réalisé dans le pays à l'époque t, nombre que nous 
avons désigné par à + &. La linéarisation habituelle donne 1, = - N'a, + p.. 
Tinbergen a trouvé, pour les U.S. A. et pour la période considérée, À' = 6, 
résultat très différent de ceux trouvés dans d'autres pays. La loi de Rueff 
donnerait (emploi inversement proportionnel au salaire réel) : 


(brad p2s8uk 
(P+p4) (a+at) 
d'où en linéarisant : 


c'est-à-dire, puisque a=1, 1=1, p= 2: 
1 
l,=-3,+ 7 Pt 
Tinbergen a lui-même trouvé (Econometrics, p. 187),( pour la Hollande 


et pour les "salaires réels" : . S , une "élasticité" par rapport à l'emploi 
t 

at qui peut varier considérablement (et même dépasser 1 en valeur abso- 

lue) suivant la politique d'exportation adoptée. 

Mais il indique qu'aux U.S.A. (peut-être en raison de la forte corréla- 
tion qui unit l'emploi a et le coût de la vie p4, les deux croissant en pério- 
de de prospérité et décroissant en période de dépression) l'influence de a, 
pour p, constant (influence traduite par le "coefficient de régression" - x 
est indécelable dans les statistiques, il prend donc : 


(10) 1, — À Pt 

il remarque d'ailleurs que la valeur extrême 0,63 qu'il trouve pour À d'a- 
près ses statistiques étant supérieure 1/2, cela signifie que la fluc- 
tuation des salaires est plus que proportionnelle à celle du coût de la vie 
(1/1 varie plus que p,/p), autrement dit, que le salaire réel (pouvoir d'a- 
chat du salaire horaire) est plus bas au cours de la dépression (période 


de baisse des prix) qu'au cours de la prospérité; cela différencie encore 
les U.S.A. de l'Europe où les salaires sont plus rigides. 


(1) Sauf dans certains marchés de gros où la spéculation sévit, comme l'a montré Four- 
geaud pour le marché du textile en France (Cahiers du Séminaire d'Econométrie, 1 - 


p. 43). 
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Vient ensuite l'équation qui détermine les profits réalisés par les en- 
treprises, les sociétés et les spéculateurs. Cette équation est fort discu- 
table. Le profit total Z + Z, par unité de temps, à l'époque t, de tous les 
entrepreneurs et spéculateurs est regardé comme l'excès du prix total 
des biens de consommation vendus (U + U,, donné par (3) et du prix total 
des biens d'investissement nouveaux produits (en écartant les biens de 
remplacement)[B +6b,,; regardé comme fluctuant dans le même sens 
que leur quantité, c'est-à-dire que le nombre d'heures de travail qui y ont 
été consacrées, 4 mois en moyenne auparavant/sur le nombre d'heures 
d'ouvriers qui ont été payées [(1 + 1,) (© + c; et b + b,,), les heures con- 
sacrées aux biens d'investissement nouveaux étant, pour la commodité 
des calculs, décomptées suivant leur montant 4 mois auparavant.] On s'é- 
tonne que les matières premières achetées à l'étranger  (quine peuvent 
certes pas être regardées comme produites par la main d'oeuvre nationale) 
ou en tout cas les fluctuations de leur valeur globale, n'intervienne pas 
dans l'équation donnée par Tinbergen. Mais il déclare (1) pp. 12 et64)que 
pendant la période considérée, les fluctuations du commerce extérieur des 
U.S. A. lui paraissent avoir eu un effet économique négligeable. Il y aurait 
là un point à élucider. 

A ces profits industriels, Tinbergen ajoute des "profits de spéculation" 
qui seraient proportionnels à la vitesse de variation des profits au cours 
de la période écoulée, soit de la forme E(z, ,-Z,,) (choisie de façon à 
simplifier les calculs.) [Ce ne peuvent être, pour l'ensemble du pays, que 
des profits nominaux et non des profits réels]. D'où, en linéarisant, le 
profit total par unité de temps à la date t (diminué de sa moyenne généra- 
lea): 

(4) 07, =00,9 p, +2 ul, et (LL, REA Zee 


(b + © étant égal à a, c'est-à-dire à 1). 


Tinbergen admet ensuite que c'est le montant des profits qui détermi- 
ne l'investissement (défini toujours par le nombre d'heures b + b, qui lui 
sont consacrées); le taux de l'intérêt à court terme m, devrait, à priori, 
avoir une influence négative sur l'investissement. Mais l'ajustement sta- 
tistique de la formule b,=&.Z,,-8.m,,a fourni à Tinbergen un coeffi- 
cient À sensiblement nul. Il remarque d'ailleurs que cela tient à ce que les 
fluctuations du taux de l'escompte, pendant la période considérée, ont été 
faibles. (Des variations beaucoup plus considérables exerceraient certaine- 
ment une influence sur les investissements). Il écrit : 


(il écarte donc le "principe d'accélération", d'après lequel les fluctuations 
de l'investissement seraient une amplification des fluctuations de la deman- 
de. Pour lui, elles reflètent surtout les fluctuations des profits). 


. Mais une partie Seulement des profits Z + 24 est distribuée aux parti- 
culiers comme revenus. Tinbergen définit E + E+ (part des profits distri- 
bués comme revenus aux non-salariés) par : 

(12) Etant 74 
La portion E; de E,qui est dépensée est regardée comme donnée par : 


(13) ES eR enr 
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ce qui suppose que la part des revenus non ouvriers, due à la spéculation, 
est dépensée après un délai plus court que le reste du revenu, majorant 
ainso E; d'une quantité proportionnelle à Z,,,- 2, (l'action de E; surles 
te ne Se manifestant d'ailleurs qu'à la date t +2, d'après la formule 


Tinbergen indique également une autre équation (14) qui rattache le 
taux de l'escompte au volume des crédits créés par les banques (crédits 
égaux à L, +E; - U,). Mais nous n'aurons pas besoin de cette équation car 
nous avons vu que, pour Tinbergen, le taux est une "variable symptomati- 
que" qui réfléchit passivement la conjoncture sans agir sur elle. 


Il est facile, entre ces 12 équations (a), d'éliminer toutes les varia- 
bles endogènes autres que l'écart de profit" Z,, et d'obtenir ainsi une 
"équation finale" qui est une équation aux différences finies ne contenant 
plus que la variable Z, et exprimant donc la valeur de Z, en fonction des 
valeurs antérieures Z,,, Z,,, etc... (Les décalages sont tous entiers en 
raison du choix fait pour l'unité de temps). (1), (3), (5) et (6) donnent, en 
effet : 


Et, =0,9pæ2u.-1. -(b, #6 


ce qui, compte tenu de (7), (9), (10) se transforme en 


(14) Es, =(059=X)r A Ten er ML 


équation qui doit être jointe à (12), (13), et à (4) que nous récrivons (comp- 
te tenu des mêmes équations que précédemment et de (11)) : 


(15) Z4=1(09æ2X) ru: tue eus +(6=1) œZ,+ E(Z: -Z,)) 


Par combinaison de (14) et de (15), u 
donnant : 


(16) Zy=Ei + DL RACINE METZ EZEE) 


Puisque (12) et (13) fournissent E; =€,€ Z. F1 (2.7 Z.), on obtient, en 
éliminant E; ,, l'équation finale : 

Zk=Ee,E2,,+0 (Ze Zi) Li CAP Li (&-1)2:-2]+8 (Zs 22) 

On résout cette équation finale (commeune équation différentielle à 
coefficients constants) en en cherchant des solutions particulières (expo- 
nentielles) de la forme Z+: = k‘, k étant une constante que l'équation finale 
assujettit à vérifier "l'équation caractéristique" : 


(17) kK-(n+a+ë)k+În+œ(i-5)+6-6,e]-0 


Les valeurs numériques des coefficients, déduites (avec fort peu de 
précision d'ailleurs) par Tinbergen de son analyse des données statisti- 
ques, le conduisent à adopter sensiblement i'équation : 


k?-1,8 k +1 —0 


qui a deux racines imaginaires k, et k,, de module 1 et d'arguments oppo- 
sés; les deux solutions particulières kt et Ke (et la solution générale 


-- et u, disparaissent tous deux 


(a) Nous avons conservé la numérotation de Tinbergen en nous passant d'une de ses équa - 
tions. 
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k! + y k£ qu'on en déduit) définissent donc des oscillations sinusoïdales 
non amorties, de période bien déterminée facilement calculable; on se 
trouverait donc exactement dans le cas limite entre l'amortissement (am- 
plitude décroissante au cours du temps) et "l'explosion" (amplitude crois- 
sante). 


Cela est moins surprenant quand on remarque que Tinbergen déclare 
avoir ajusté ses coefficients de façon à arriver à ce résultat (qui lui paraît 
vraisemblable au vu de la sensibilité de l'économie américaine pendant la 
période considérée). 


Il est intéressant de remarquer que, dans une autre publication, un 
calcul analogue conduit sur l'économie anglaise de 1870 à 1914 lui a donné 
l'équation caractéristique : 


k?-1,27k +0,6—0 


dont les racines, cette fois toutes deux imaginaires de module < 1, condui- 
sent à un mouvement amorti. 


Tinbergen se livre ensuite à un calcul de la période des oscillations de 
Z. et de celles des oscillations supplémentaires que peuvent présenter les 
variables u,, 1,, p., etc... Ce calcul présente peu d'intérêt, étant donné 
que les "périodes propres" du système économique ne seront presque ja- 
mais observables, étant détruites par les influences exogènes et les per- 
turbations aléatoires. 


Mais le point intéressant de l'étude porte sur les mesures de "politi- 
que économique" que coit prendre l'Etat pour modifier les coefficients de 
l'équation de façon à augmenter l'amortissement, c'est-à-dire à stabiliser 
l'économie. La fixation des prix (action sur le coefficient x) et des salai- 
res (action sur À ) apparaissent, dans le modèle établi par Tinbergen, 
comme sans influence, puisque les coefficients À et rt n'apparaissent pas 
dans l'équation finale (a). 


Très grande, au contraire, apparaît l'influence de la spéculation qui 
se manifeste par les coefficients n et Ë, ainsi que celle (traduite par le 
coefficient «) de la sensibilité des investissements par rapport aux profits. 
L'amortissement sera d'autant plus grand que les racines de l'équation (17) 
seront plus petites en module, c'est-à-dire que la somme k, + k, et le pro- 
duit k,k, des racines seront plus petits en valeur absolue. Il y a donc inté- 
rêt à rendre les plus petits possibles, dans (17), le coefficient de k et le 
terme constant, ce qui sera réalisé en rendant les plus petits possibles les 
coefficients n et E (ce qui revient à freiner la spéculation par la réglemen- 
tation du crédit ou par un contrôle boursier), les coefficients €, et € (ce 
qui revient à freiner les fluctuations E; des revenus non-ouvriers dépen- 
sés (b) et le coefficient &« (qui peut être maintenu nul en réalisant la cons- 
tance des investissements à l'aide d'investissements compensateurs effec- 
tués par l'Etat). 


(a) Cela tient à ce que les salariés sont à la fois producteurs et consommateurs.Toutefois 
le niveau des prix et des salaires, s'il est combiné avec une dévaluation, a une influ - 
ence sur l'exportation et l'importation, dont Tinbergen a fait abstraction pour les fluc- 
tuations aux U.S. À., mais qu'il a étudié en détail pour les Pays-Bas. 


(b) Mais l'étude de l'incidence de la, politique fiscale devrait comme H. Eyraud enparticuw 
lier l'a montré(3) prendre en considération sa profonde répercussion sur la production 
en remaniant les équations utilisées. 
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. Tinbergen fait d'ailleurs remarquer qu'une politique qui ne s'attaque- 
rait pas à la spéculation et qui serait centrée sur la seule stabilisation des 
investissements, c'est-à-dire qui aurait pour objectif unique de réaliser 
%= 0, serait dangereuse. En effet, le coefficient à (qui traduit le profit 
retiré de la production de biens d'investissement) a dans les statistiques 
de Tinbergen une valeur numérique qui, quoique assez incertaine, est net- 
tement > 1. Il en résulte q'une diminution de & se traduit par une augmen- 
tation du produit des racines de (17) qui peut alors, si la spéculation n'est 
pas jugulée, dépasser la valeur 1, ce qui fait tomber dans le cas d'instabi- 
lité, dans le cas "explosif". 


Toutes ces remarques sont des exemples des nombreux problèmes que 
l'on peut résoudre sur un modèle économique déterminé. Il s'agit ici d'un 
problème ultra-simplifié, en particulier du fait que l'on n'a pas pris en 
considération le commerce extérieur et le change, ce que Tinbergen juge 
inutile pour les U.S.A. Il a d'ailleurs remanié son modèle pour l'adapter 
à d'autres pays comme la Hollande où le commerce extérieur est essen- 
tiel; la discussion est alors plus complexe en ce qui concerne l'incidence 
de la politique des salaires et des prix (c'est pourquoi je la laisse de côté 
pour ne pas alourdir l'exposé) quoiqu'elle n'introduise aucune difficulté 
mathématique nouvelle. Ce que je veux souligner pour terminer, c'est que, 
si la résolution d'un modèle mathématique donné ne présente aucune diffi- 
culté de principe (en tout cas en ce qui concerne les petits mouvements 
pour lesquels l'approximation linéaire est valable), l'obstable redoutable 
auquel se heurtera longtemps l'économétrie réside dans l'insuffisance des 
preuves statistiques contraignantes en faveur detelle outelle des équations 
de structure proposées par les économistes, dans l'importance des pertur- 
bations aléatoires qui modifient les variables, et dans la grande impréci- 
sion avec laquelle sont connus les coefficients des équations (sans parler 
des valeurs d'équilibre et des tendances séculaires dont la détermination 
pose aussi des problèmes délicats). L'amélioration et la normalisation du 
matériel statistique, sur lesquelles l'attention des gouvernements commen- 
ce à se porter, sont une oeuvre de longue échéance; c'est pourtant la base 
indispensable pour la construction d'une "économie mathématique" qui soit 
à l'épreuve des faits. 
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